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1 Introduction 

The classical Kirchhoff Approximation (KA) in diffraction theory is an asymp- 
totic relation between the Dirichlet and Neumann data of a solution of an 
exterior boundary value problem for the Helmholtz equation with frequency 
parameter k ^> 1. The KA is sensitive to the assumption of positive curvature 
of the boundary and its accuracy deteriorates in the presence of flattening 
regions |J|. Numerical methods for high-frequency problems have attracted 
much attention lately - see e.g. j2] • In an attempt to include small and vanish- 
ing curvatures uniformly in an asymptotic theory, we propose to study high- 
frequency asymptotic properties of the Dirichlet-to-Neumann (DtN) operator 
rather than those of an individual solution. In j^j, we conjectured that the 
pseudodifferential symbol of the DtN operator, appropriately scaled, tends 
to a simple universal function as k —>■ oo. Here we report results of a more 
detailed numerical study. The results support the said conjecture in the case 
of a convex scatterer. In a non-convex case, we observe a deviation from the 
universal limit function in a narrow range of it's argument values. 
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2 The DtN operator and the Limit Shape 
Hypothesis 

Consider the Helmholtz equation Au + k 2 u = in the exterior of a simply 
connected bounded domain Q C R 2 with smooth boundary T. Given a 
function / on T of a certain regularity [HUE], the Dirichlet problem u\r = 
f has a unique solution u satisfying the Sommerfeld radiation condition 
d r u — iku = o(r -1 / 2 ) as r — > oo. The normal derivative g = d u u\r is 
a function of known regularity on T. The map Af : / — > g is called the 
Dirichlet-to-Neumann (DtN) operator. 

Let s be the arclength parameter on T, and L the length of T. Set 
<p = 2tts/L. The Dirichlet and Neumann data / and g in the above diffraction 
problem are 27r-periodic functions of variable <p. Let f(4>) = ^f{n)e m ^ be 
the Fourier series of /. Write M as a periodic pseudodifferential operator 
(PPDO) [UE1 

oo 

g(<f>)=Mf(<f>) = <7(<f>,n)f(n)e in <t>. (1) 

n=— oo 

The function a(4>, n) = e - %n ^ ]\[ e m <i> j s called the symbol of N '. 

The operator M depends on the boundary V as well as on the frequency 
k. We reflect this in notation of the symbol by writing o~(<p, n) = <7r(0, n; k). 
In j^j we noted an universal (r-independent) high-frequency asymptotic be- 

haviour of the symbol as a function of variable £ = £(n, k) = — — . Define 



Lk 



O"lim(0 




KI<1, 
KI>1. 



Hypothesis 1. For any e > and am/ > 1, i/iere exists k* > suc/i t/iat 

AT 1 ot(0, fe) - cum (£(n, fc)) <e (2) 

whenever k > k* and \£(n,k)\ < £*. 

Here are some theoretical arguments in favour of Hypothesis 1. 

1. The statement holds if T is a circle of any radius jHj- 

2. If |£(n, fc)| > 1, then the inequality (J2J can be established by construct- 
ing an asymptotic WKB solution, as pointed out by L. Friedlander 
(Univ. of Arizona), personal communication, February 2004. 
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3. A simple if not completely rigorous argument shows that the hypothesis 
is consistent with KA for a convex domain [5]. 

Yet we admit that Hypothesis 1 may be true for some classes of boundary 
curves and false for others. To restore status quo with numerical experiment, 
we formulate a somewhat weaker Hypothesis 2 below. 

Note that the symbol ar((p,n; k) generally depends on (except when 
T is a circle), while the limit function is 0-independent. So we are trying 
to approximate the DtN operator by a shift-invariant PPDO. It can only be 
possible if the Fourier series of the symbol in asymptotically reduces to a 
single constant term. Put (omitting the subscript V in the right-hand side) 

a T ((f>, n- k) = a (n; k) + a ±1 (n; k)e ±i(t> + a± 2 (n; k)e ±2i(t> + ■■■. (3) 

We shall compare the mean symbol cr$(n] k) = (2n)^ 1 J 27r 0"r(0, n; k) d<fi to 
the limit function and watch whether the / 2 -norm ||<3"'(n; k)\ \ of a bi-infinite 
vector formed by the rest of Fourier coefficients Q is relatively small. Recall: 

\\a'(n;k)\\ 2 =Y,\t m (n;k)\ 2 . (4) 

Hypothesis 2. For any boundary curve T and any given £* > 1, e > 0, and 

5 > 0, there exists k* > such that if k > and \C,(n, k)\ < £*, then 

1. the shape of the mean symbol cr of the DtN operator follows that of 

0"lim-' 

\k~ l a (£(n, k); k) - a hm (0 | < e; (5) 

2. the remaining Fourier coefficients of the symbol are collectively small: 

k~ l \\a'{n, k)\\ < e, (6) 

if distance (£(n; k), I) > 5. Here I is either the empty set or a certain 
"exceptional" set determined by the curve T. 

Note that Hypothesis 1 implies Hypothesis 2 with 7 = 0. The parameter 
5 in Hypothesis 2 is introduced to account for a non-uniform convergence 
near I when I is nonempty. Note also that in this paper we require T to be 
a smooth curve, but there exist numerical results supporting validity of the 
statement for domains with corners. 
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3 Methodology of numerical verification 



To test the hypothesis numerically, we use known sample solutions satisfy- 
ing the Helmholtz equation (HE) in the exterior domain R 2 \ Q and the 
radiation condition (RC), and compute Fourier coefficients of the Dirichlet 
and Neumann data. Solutions of HE in R 2 \ (0, 0) with wavenumber k and 
satisfying RC are spanned by the Hankel functions H$(kr), m — 0, 1, • ■ •, 
r = \f\ = yjx' 1 + y 2 . The origin can be viewed as an emitter, or source. 
Now, by taking fictitious sources at arbitrary locations S G Q, the family 
H$(k\r — S\) of sample solutions in R 2 \ f2 is constructed. For the verifica- 
tion procedure one can use a countable sub-family with linear combinations 
dense in the space of solutions. In this work, we use i/p -solutions with 
sources near the boundary T and approximately equidistributed along T. A 
possibility to represent an arbitrary solution of HE+RC in the form of a 
single layer potential (provided k 2 is not an interior eigenvalue |3 §3.2.1]), 
justifies this choice. An extreme opposite possibility is to choose a family 
of if^-solutions, m = 0, 1, 2, • ■ -, with fixed source. It needs the Rayleigh 
hypothesis for domain Q to hold, which is true, for example, if T is an ellipse 
with eccentricity e < 1/V2 [TH] . 

Let us first describe a procedure used in 0. Take a uniform partition 
{Pi}, I — 1, 2, ... , Z max of the curve T. Evaluate a sample solution Ho~\k\r — 
S\) and its normal derivative at the points f = Pi to obtain the vectors fi and 
gi of size / max . Then compute the discrete Fourier transforms and consider 
their truncations f(n), gin), \n\ < n max . Find the ratio a{n) = f{n)/g{n) 
and compare k~ l <r(n) to o"i im f^rj to verify Hypothesis 1. Typically in our 
examples Lk ~ 10 2 -r- 10 3 ; we chose Z max = 2 12 -=-2 24 , and n max w 3Lk. Higher 
Fourier coefficients are vanishingly small, that is why we cut them off. 

In more detail, let r (</>), <p [0, 27r] be the parametrization of T by the 
normalized arclength = 2ns /L. Put (pi = 27r/// max and Pi = r(4>i). Then 

/, = H^ikm) - S\), gi = -ei k [l - (r'(0O) 2 ] V2 H^iklf^i) - S\). 

Here r(0) = \r{<p)\ and Ei = (— 1) M , where pi is the number of intersections 
of T with the open interval {SPi). Note that pi = and E\ = 1 if Q is convex. 

Our judgement about validity of Hypothesis 1 in [3] was based on the 
outlined procedure, where we effectively kept over the mean symbol <5"o(n; k) 
only. But this is not enough. Let us engage in the study of components of 
the vector a', see (HJ), - apart from <7 (n; k). Now we take several sources, 
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Si, . . . , Sj at once. Assume, for the sake of symmetry, that J is odd, J = 
2w m ax + 1- Denote by // and gj the data of the solution with source at 
Sj, and by f^{n), g^(n) the corresponding components of the (truncated) 
discrete Fourier transforms. The following relations follow from Q,®: for 
every j = 1, . . . , J 

g>(n) = ^2 f (n - m) a m (n - m; k). (7) 

m 

Reduce the infinite summation to a finite number of terms keeping only the 
components a m ( ■ ; k) with |m| < m max . For example, if J = 3, then for each 
n = — Timax, • • • , «max after cut-off we get a linear system of three equations 
with three unknowns o"o(n; k), <r±i(n; k) : 

gi(n) = f j (n + l)a_ x {n + 1) + f (n)a (n) + f j (n - l)a x {n - 1), j = 1, 2, 3. 

Solving all obtained systems, we approximately find a m {n) for (at least) \n\ < 
n max — m max and \m\ < m max . Now the left-hand sides of the inequalities (jHJ), 
© can be evaluated; of course, summation in (jlj) is restricted to \m\ < m max . 




(a) -0.5 0.5 (b) -10 1 



Figure 1: Test domains: (a) Convex (ellipse), (b) Non-convex (kite jl]) 



4 Results and discussion 

We present results for two symmetric domains shown on Fig. 1: the el- 
lipse x(t) = cost, y(t) = 0.6 sin t, and a non-convex "kite" x(t) = cost + 
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0.65 cos2t — 0.65, y(t) = 1.5 sint. If = at the right x-intercept of T, then 
due to symmetry a m (—n; k) = a_ m {n, k) and it suffices to study the symbols 
for n > 0. 

The real and imaginary parts of the rescaled mean symbol k~ 1 ao{n, k) 
are compared to the limit curves on Fig. 2, 3. Here £ = k) as defined in 
Sect. 2. The parameters are: frequency k = 200; number of sources J = 201. 
The kite's curves exhibit some roughness when £ e (0.8, 1). 




(a) 0.5 1 1.5 ' (b) 0.5 1 1.5 



Figure 2: — Re<T (n; k)/k and — Reo"i im (£) vs £ = k) 




200 



(a) 0.5 1 1.5 (b) 0.5 1 1.5 

Figure 3: lma (n;k)/k and Im<7i im (£) vs £ = 

Fig. 4 shows the left-hand side of the inequality © vs £(n, k) for /c = 50, 
100, 200, 400, 800. The value of J was always set equal to k + 1. In the 
case of ellipse, the norm shrinks to naught as k grows. It isn't quite so for 
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Figure 4: Violation of symbol shift-invariance: k 1 \\a'(n; k)\\ vs £(n, k) 



the kite. The peak over the interval (0.8, 1) stays steady. In the frameworks 
of Hypothesis 2, we say that the exceptional set / is empty for the ellipse, 
though the convergence near |£| = 1 is much slower than away from |£| = 1. 
The set I for the kite is apparently contained in the union (—1, — 0.8)U(0.8, 1). 

Computational note. Computation of the Fourier coefficients <x m (n; k) 
of the symbol requires solution of truncated systems (|7jl . If the cutoff sub- 
script is rather large, one has to take trouble to ensure that Fourier coeffi- 
cients P(n — m) are not vanishingly small. To this end, the sources should 
be placed close to the boundary, preventing the Dirichlet data of sample so- 
lutions from being "too smooth". The reported results are obtained with 
sources located at the distance from about 10~ 2 to 10~ 3 (for larger values of 
k) from T. 
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